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Abstract 

Although we have yet to determine whether the DAMA data represents a true 
discovery of new physics, among such interpretations inelastic dark matter (IDM) can 
match the energy spectrum of DAMA very well while not contradicting the results 
of other direct detection searches. In this paper we investigate the general properties 
that a viable IDM candidate must have and search for simple models that realize these 
properties in natural ways. We begin by determining the regions of IDM parameter 
space that are allowed by direct detection searches including DAMA, paying special 
attention to larger IDM masses. We observe that an inelastic dark matter candidate 
with electroweak interactions can naturally satisfy observational constraints while si- 
multaneously yielding the correct thermal relic abundance. We comment on several 
other proposed dark matter explanations for the DAMA signal and demonstrate that 
one of the proposed alternatives - elastic scattering of dark matter off electrons - is 
strongly disfavored when the modulated and unmodulated DAMA spectral data are 
taken into account. We then outline the general essential features of IDM models in 
which inelastic scattering off nuclei is mediated by the exchange of a massive gauge 
boson, and construct natural models in the context of a warped extra dimension and 
supersymmetry. 
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1 Introduction 



The DAMA/Nal and DAMA/LIBRA experiments observe an annual modulation signal in 
their Nal-based scintillation detectors with a statistical significance of 8.3 a [Ij. A possible 
origin of this signal is galactic dark matter (DM) scattering off the nuclei in the detectors 
of these experiments [2|,[3]. The annual variation would then result from the motion of the 
Earth relative to the Sun as it passes through the halo of dark matter enveloping our galaxy. 
The phase, period, and amplitude of the modulation signal seen by DAMA are all consistent 
with DM scattering [Ij. We have yet to determine whether the observation is the result of 
truly new physics, or if it results from an unaccounted-for detector effect or background. 
In this paper, we ask the question: what if the DAMA result is truly a discovery of dark 
matter? What could it possibly be? 

The major challenge for a DM interpretation of the DAMA result is that it appears 
to be at odds with the bounds on coherent DM-nucleus scattering obtained by other DM 
direct detection experiments such as CDMS [1] and XENON [5J. In contrast to DAMA, 
these experiments search for an unmodulated nuclear recoil signal from DM scattering using 
stronger background rejection methods. The bounds they place on DM-nucleon cross sections 
rule out coherent elastic scattering off iodine nuclei as the origin of the DAMA signal by 
several orders of magnitude. Lighter DM scattering coherently off sodium nuclei is marginally 
consistent with both DAMA and other experiments [6l [7], but gives a very poor fit to the 
energy spectrum of the modulated DAMA signal [H [9l [TO] . 

An elegant possibility that can account for the signal observed by the DAMA experiments 
that is also consistent with other direct detection experiments is dark matter that scatters 
inelastically off nuclei [TTl [T2| [T3| [T^. In the inelastic dark matter (IDM) scenario, the dark 
matter particle Xi scatters preferentially off target nuclei into a slightly heavier X2 state. 
The kinematics of this process can enhance the nuclear recoil signal at DAMA relative to 
other experiments such as CDMS in a couple of ways. To produce a nuclear recoil signal 
with energy Eji, the minimum incident velocity of the DM particle is [TT] 



where 6 is the mass splitting between xi and X2, is the mass of the target nucleus, 
and /iTv is the reduced mass of the nucleus-DM system. The distribution of DM velocities 
in the galactic halo is expected to be approximately Maxwellian with an upper cutoff at 
the galactic escape velocity Vesc [3]. When the second term in this expression dominates, 
the minimal velocity needed to produce a recoil energy is lower for heavier nuclei. This 
leads to an enhanced signal at DAMA, which contains iodine with A ~ 127 as a detector 
material, relative to CDMS, consisting of germanium with A ~ 73 [llj. The kinematics of 
inelastic DM scattering also increases the amount of annual modulation compared to the 
unmodulated signal rate, further enhancing the signal at DAMA relative to other direct 
detection experiments [llj. 

While IDM provides a compelling explanation for the DAMA signal, only a few concrete 
particle physics candidates have been proposed. A model with sneutrino DM where the 
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inelastic splitting is induced by the lepton-number violating superpotential operator W D 
{L- / A was suggested in Ref. [11]. Ref. [15] proposed that the inelasticity could arise 
from the radiative splitting of masses within a multiplet after the spontaneous breakdown 
of a new non-Abelian hidden gauge symmetry around a GeV. Pseudo-Dirac neutralinos 
as IDM arising from from approximately i?-symmetric SUSY scenarios are considered in 
Refs. [Ml HI]. 

In the present work we seek to obtain a broader overview of potential candidates for 
IDM. We study the general features required for IDM to account for the DAMA signal and 
we describe several explicit IDM candidates. To remain as general as possible, we do not 
attempt to account for the tantalizing indirect hints for dark matter such as the excess 
positron and electron fluxes observed by PAMELA [18], ATIC [19], and PPB-BETS [20j, 
the INTEGRAL 511 keV line [21], or the WMAP haze [221 [231 [2l] . In this paper our goal 
is simply to understand the possibilities for DAMA alone. What if DAMA represents a true 
discovery? What would be credible candidates for what it could be? It would of course be 
interesting to determine which of the candidates for IDM might be compatible with these 
indirect signals, but we postpone this direction to future work. 

The outline of this paper is as follows. In Section [2] we investigate how well the hypothesis 
of IDM can account for the DAMA signal while evading the constraints from other direct 
detection searches for dark matter. We also argue that several other mechanisms proposed to 
explain this puzzle do not give a good fit to the full DAMA dataset. In Section [3] we discuss 
the general features required for a model to generate IDM, and we study the corresponding 
phenomenological constraints on these features. We present several plausible and explicit 
models that can give rise to IDM in Section |H Section [5] is reserved for our conclusions. 

2 Inelastic Dark Matter as an Explanation for DAMA 

Inelastic dark matter has been shown to provide a compelling explanation for the signal in 
the DAMA experiments while remaining compatible with other direct DM searches such as 
CDMS, XENON, and CRESST for dark matter masses at least as large as 250 GeV [H]. 
Here we extend the analysis of Ref. ^Mj to larger dark matter masses. Such larger masses 
are of interest theoretically because they can be natural in the context of models addressing 
the hierarchy problem such as warped geometry p5j and supersymmetry [26]. In the case 
of electroweakly charged dark matter, they also lead to the correct thermal relic abundance. 
Additional attention to heavy dark matter is due to the recent experimental results from 
ATIC and PPB-BETS. We also consider some of the proposed alternative explanations of 
the DAMA signal and comment on their viability. 

2.1 IDM Fits to the DAMA Data 

We begin by reviewing the formalism for calculating the expected signal at direct detection 
experiments from the scattering of dark matter. The total rate of inelastic nuclear recoil 
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scatterings per unit mass of detector per unit recoil energy Eji in the lab frame is 

-Nr-^ f (2.i; 



where pdm — 0.3 GeV/cm^ is the local DM density, M^m is the DM mass, is the number 
of target nuclei per unit mass of detector, and f{v,Ve) is the local dark matter velocity 
distribution. For coherent spin-independent DM scattering, the DM-nucleus differential 
cross section da/dE^ has the general form [271 [28] 

da _ I niN al [fp Z + fn{A- Z)f . 



AT? 2 O 2 f2 ^ ^^r)^ (2-2) 

where mjv is the mass of the target nucleus with atomic and molecular numbers A and Z, 
[In is the DM-nucleon reduced mass, /p and are effective coherent couplings to the proton 
and neutron, and cr^ is the overall effective DM-neutron cross section at zero momentum 
transfer. The function F'^{Epi) is a form factor characterizing the loss of coherence as the 
momentum transfer = 2m]sf Eji deviates from zero. For computational simplicity, we use 
the Helm/Lewin-Smith [291 EQ] parameterization of the form factor. In order to correct for 
the fact that this parameterization can be off by ~ 20% for larger values of Eji, we weight 
this by a quartic polynomial fit to the table given in [31], which gives the ratios of the 
Helm/Lewin-Smith form factor to the more accurate Two- Parameter Fermi (Woods-Saxon) 
form factor for various elements and values of Eft. 

We take the DM velocity distribution to be Maxwellian with a cutoff [3], 

= r^Wl^"^''^'^^^'^'" e{v,sc-\v + Ve\). (2.3) 

Here, Vq ~ 220 km /s is the DM rms speed, Vesc is the local DM escape velocity in the halo 
frame, Ve is the velocity of the Earth with respect to the galactic DM halo, and v is the DM 
velocity in the Earth frame. The signals from IDM are extremely sensitive to the value of 
Vesc, which is thought to be in the range 498 km/ s < v^sc < 608 km/ s [32]. The velocity of 
the Earth relative to the halo, Vg, has components from both the motion of the solar system 
relative to the halo as well as the annual motion of the Earth about the Sun, 

Ve{t) =v^ + Vo [ei cos(27r(t - ti)) + sin(27r(t - t^))\ (2.4) 

where ~ (0,220,0) + (10,5,7) km/s is the velocity of the solar system relative to 
the halo [33], |3lj, Vo — 29.79 km/s is the Earth's orbital speed [30], and t is measured 
in years. Following the conventions of Ref. [35] and the discussion in Ref. [10], we are 
using coordinates where x points to the center of the galaxy, y gives the direction of disk 
rotation, and z points to the north galactic pole. The directions of the Earth's motion 
on ti = March 21 (ei) and June 21 (eg) are given by ei = (0.9931,0.1170,-0.01032) and 
62 = (-0.0670,0.4927,-0.8678) [35]. 

The annual variation of the scattering rate due to Ve{t) is very nearly sinusoidal and we 
estimate the amplitude of the modulated rate to be 
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When computing the unmodulated rate for a given experiment, we integrate the total rate 
over the time periods that the experiment recorded data. To convert rates to detector 
signals, we rescale by the efficiency of the detector and account for quenching, but we do 
not include any detector resolution effects. In addition, we have included a correction for 
channeling effects at DAMA, but find them to be unimportant for the range of parameters 
we consider^ 

To determine the extent to which IDM can account for the DAMA signal, we compute 
the modulated recoil spectrum of IDM candidates and compare this to the energy spectrum 
in the twelve lowest bins in the 2—8 keVee range reported by DAMA in Ref . [1] . For a given 
dark matter particle mass, we use a goodness-of-fit metric to determine the 90% and 99% 
confidence level allowed regions. We define this metric as 



where denotes the average of the left, center, and right values of S in the i*'* 0.5 keV 
width bin, and cr^^ata reported uncertainty in the measurement. For each value of the 
dark matter mass, we scan over 2 parameters (the overall nucleon cross section ct°, and the 
mass splitting 5), so we require that < 16.0 (23.2) at the 90 (99)% level for 12 - 2 = 10 
degrees of freedom. We only consider signals from the scattering off of iodine, as it is expected 
to completely dominate for the parameters of interest to us. Following Ref. [H], we take the 
quenching factor for iodine to be qi = 0.085. 

As well as fitting to the DAMA results, we also compute the signals that each IDM can- 
didate would produce at the CDMS experiments [371 [Mill], CRESST-II [39], and ZEPLIN- 
III [10]. This imposes further constraints on the properties of a potential IDM candidate|§ 
In order to be conservative in excluding IDM parameter space, we will assume that the small 
number of events seen by these experiments are signal events, and use Poisson statistics to 
find the region of parameter space excluded at the 99% confidence level based on the number 
of observed events. 

CDMS has published data from three runs at the Soudan Underground Laboratory with 
approximate exposures of 19.4 kg-day [37], 34 kg-day ^38j, and 121.3 kg-day [3]. In total, 
these experiments reported two events between 10 keV and 100 keV, which we assume to be 
signal. We thus require that the expected total number of events, integrated over the time the 
experiments ran, obeys Ntot < 8.4 at the 99% confidence level. We only consider scattering 
off germanium as it is expected to dominate for the heavier dark matter we consider. 

CRESST-II has published data from a run in 2004 using prototype detector modules [13] 
and more recently has published results from a commissioning run carried out in 2007 |39j . 
Since there were significant changes to the detector modules between these runs, including 

^The fraction of channeled events falls quickly with recoil energy, and can be approximated for iodine as 
// ~ io~V^W(ii^5kcV) [2gl[TUj. Since IDM suppresses scattering processes with low Eji, one expects only a 
very small number of channeled events. 

^As in the analysis of Ref. [M], we find that the constraints from other experiments such as XENONIO [5], 
KIMS [H], and ZEPLIN-II [42] are currently not as important so we have not included them in our plots. 




(2.6) 
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the addition of neutron shielding, we opt not to combine these data sets and consider 
constraints only from the later commissioning run, which had an exposure of 47.9 kg-day 
and an acceptance of 0.9 for tungsten recoilsjj Taking the seven observed events between 
12keV and 100 keV to be signal, we require that the total number of predicted tungsten 
recoil events, integrated over the duration of the experiment, obeys Ntot < 16.0 at the 99% 
confidence level. 

The ZEPLIN-III experiment has recently released data from a run in 2008 with an 
effective exposure of 126.7 kg-day [40j. The experiment observed seven events in its liquid 
xenon detector between 2 keVee and 16 keVee, so we require that the total number of events 
obeys Ntot < 16.0 at the 99% confidence level. To convert measured energy to recoil 
energy E^, we use the energy- dependent quenching factor given in Fig. 15 of Ref. [ID], which 
saturates at qxe — 0.48 around E^ = lOkeVee. Below this scale, we use the parametrization 
qxe ~ (0.142 Ed + 0.005) Exp [-0.305 E/^^^] given in Eq. (4.3) of Ref. [g, which we find 
gives a good fit to the curve. 

In Figs. [1] and [2] we show the region allowed by DAMA and constraints from CDMS II, 
CRESST-II, and ZEPLIN-III at the fiducial point = fp for various values of the dark 
matter mass between 100 GeV and 5 TeV. (Note that ultimately the relation between /„ 
and fp is model-dependent, but we make this choice to allow a straightforward comparison 
to previous work [6l[7l[TTl[T2l[T3l[Tlj.) In addition, we consider escape velocities of 500 km/s 
and 600 km/s. We observe that while the relative constraints from CDMS II, CRESST-II, 
and ZEPLIN-III become stronger for heavier dark matter candidates, very heavy dark matter 
is not ruled out. While both the required cross section for DAMA and the constraint curves 
move upwards for heavier dark matter, they do so at roughly the same rate. One can 
understand this asymptotic behavior as follows. For large values of the dark matter mass 
relative to the mass of the nucleus, the reduced mass in Eq. (11.11) becomes /i^r ~ m^r and is 
independent of the dark matter mass. Therefore the only mass dependence is in the prefactor 
of Eq. (12. ip which is the same for all experiments, so the constraint curves from different 
experiments do not move relative to each other. 

It is quite interesting that heavier dark matter may be allowed - this opens up the 
possibility of simple dark matter models for which heavier masses are preferred for getting 
the right thermal relic abundance, or for explaining indirect signals such as the e"*" + e~ 
excess at ~ 300 — 800 GeV seen by ATIC and PPB-BETS. Models addressing the hierarchy 
problem also often prefer heavier dark matter candidates [251 126] . 

We stress, however, that the results of the present section depend on a number of 
astrophysical and nuclear physics quantities that are not fully understood and have large 
uncertainties. The size of the modulated signal at DAMA is particularly sensitive to the local 
dark matter velocity distribution. For example, as one can see from Figs. [1] and [21 increasing 
the halo escape velocity from 500 km/s to 600 km/s significantly tightens the constraints 
from CDMS. This is because inelastic DM scattering requires one to sample from the tail of 
the velocity distribution, especially for lighter target elements. Deviations from the assumed 

■^We note that the combined data sets would allow parameter points that are ruled out according to the 
later commissioning run considered in isolation due to the large number of observed events in the earlier run 
(5) relative to its exposure (20.5 kg-day). 
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Figure 1: Allowed regions for fits of inelastic dark matter to the DAMA data, as well as 
constraints from CDMS II, CRESST-II, and ZEPLIN-III. We fix the local DM density at 
p = 0.3GeV/cm^, and vary the DM mass and escape velocity. These plots assume the 
relation fp = fn- 
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Figure 2: Allowed regions for fits of inelastic dark matter to the DAMA data, as well as 
constraints from CDMS II, CRESST-II, and ZEPLIN-IIL We fix the local DM density at 
p = 0.3GeV/cm^, and vary the DM mass and escape velocity. These plots assume the 
relation fp = fn- 
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Maxwellian distribution at high velocities, as well as additional DM substructures such as 
streams or sub-halos, can also significantly affect the allowed region [l5l HBl HZl Hi]. For 
this reason, we have been extremely conservative in identifying the excluded region by using 
a goodness-of-fit estimator and only showing 99% confidence level exclusion contours for 
CDMS II, CRESST-II, and ZEPLIN-III. Given the uncertainties, it would be premature to 
rule out a broader range of IDM parameter space. 

A second important uncertainty is the value of the local dark matter density pdm, which 
is known only up to a factor of ~ 2. Varying p_Dj\/ will affect the overall normalization of the 
cross section 0"° in the plots of Figs. [1] and [2] by a factor inversely proportional to it. While 
this is not terribly important for demonstrating the existence of an allowed region, since it 
affects all signals equally, the uncertainty in p^A/ can be very important when comparing 
the allowed region to a model that predicts a specific value for cr°. 

Finally, we wish to emphasize that the analysis performed in this section (and in Refs. 
[inilll]) is not model-independent. Simple models of IDM, for example dark matter charged 
under SU{2)l, will not respect the relation fp = fn, and the curves in Figs. [1] and [2] will move 
around by different amounts that depend on the atomic numbers relevant to the experiment. 
In Section [3l we will redo the analysis for dark matter charged under SU{2)l, which scatters 
through Z°-exchange. There we will see that when one considers the heavier dark matter 
masses that yield the correct thermal relic abundance [H], the scattering cross section from 
Z°-exchange is tantalizingly close to the DAMA preferred region (which remains qualitatively 
similar to the fp = fn case considered here). 

2.2 Other DM Explanations for DAMA 

In addition to IDM, a number of other non-standard DM candidates have been suggested as 
possible sources for the DAMA signal. These include lighter spin-independent elastic DM [6j, 
lighter spin- dependent elastic DM [59], and elastic scattering off atomic electrons ^50j. These 
alternatives to IDM turn out to be extremely strongly constrained when spectral data from 
DAMA is included in the analysis p[9l[Tn]. 

Light 2 — 10 GeV elastic, spin-independent DM was proposed as an explanation for the 
DAMA signal in Ref. [6]. The DAMA signal in this scenario comes primarily from the light 
DM scattering off sodium rather than iodine. Such a light DM state would produce recoils 
near the lower end of the sensitivity range of most other direct detection experiments which 
consist of heavier target nuclei, strongly suppressing their signals relative to DAMA. The 
effect of light DM on DAMA can be further enhanced relative to other direct detection 
probes by channeling, which effectively reduces the amount of quenching in the DAMA 
target, leading to more events in the range of sensitivity [7]. 

However, subsequent analyses taking into account the modulated and unmodulated 
single-hit DAMA energy spectra indicate that this light DM scenario is strongly disfavored. 
It is found that light DM either does not provide a good fit to the modulated spectrum or 
predicts an unmodulated single-hit rate in the lowest energy bins that is much larger than 
the total rate observed by DAMA [8]. Even so, light spin- dependent DM may still be viable 
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if channeling is included [TO]. 

A second alternative explanation for the DAMA result consists of moderately heavy 
{rriDM ^ 10 GeV) elastic DM that scatters primarily off atomic electrons rather than 
nuclei [50]. The vast majority of these scatterings would produce an electromagnetic signal 
in the DAMA detector in the eV energy range, well below the keV energies to which this 
detector is sensitive. However, the scattering of a halo DM particle off an atomic electron 
with an unusually high momentum in the tail of its distribution, on the order of an MeV, can 
generate a detected electromagnetic signal at DAMA of -E^ ~ few keV [50]. No such signal 
would have been recorded in other direct detection experiments such as CDMS and XENON 
since these experiments are careful to filter out electromagnetic events that they expect to 
arise from backgrounds. Electron scattering DM is also attractive in light of the PAMELA 
results [STl [52] , which can be interpreted as coming from DM annihilating preferentially into 
leptons [531 Eg. 

Following the analysis of Ref. [8j for light DM scattering off nuclei, we investigate whether 
electron interacting DM is consistent with the modulated and unmodulated single-hit DAMA 
energy spectra. In our analysis we compute the modulated and unmodulated rates for DM 
scattering off electrons as in Ref. ^50j. We compare our binned results to the lowest twelve 
2-8 keVee DAMA modulated bins and the lowest six 0.875-2.125 keVee DAMA unmodulated 
single-hit bins using a modified measure. We use a standard goodness- of- fit measure for 
the modulated bins, while for the unmodulated bins we add to the only if the predicted 
signal is larger than the observed value to allow for an unmodulated background. This 
procedure is very conservative in that it will only underestimate the excluded regions. 

We find that under the assumptions about the DM made in Ref. [50], namely that the 
DM is fermionic and interacts with quarks by the exchange of a scalar or a gauge boson 
with {V ± A) couplings]^ electron-interacting DM as an explanation for DAMA is excluded 
well beyond the 99% confidence level. This occurs for precisely the same reasons that light 
elastic DM is strongly disfavored: either the modulated signal is too low, or the unmodulated 
single-hit signal {i.e., the signal excluding multiple scintillation events) exceeds the total rate 
observed by DAMA. In the present case, the signal rate falls quickly with increasing detected 
energy Ed because the momentum distribution of atomic electrons decreases rapidly in the 
relevant range, approximately as |50j . 

This tension is illustrated in Fig. [3] where we show the best fit (lowest effective x^) to 
the full DAMA spectral dataset (Fit A), as well as the best fit (lowest effective x^) to the 
modulated dataset alone leaving out the lowest 2 keVee energy bin (Fit B). In making these 
fits, we assume either chiral vector {V ± A) or scalar four-fermion interactions between the 
electron and a fermionic DM particle of mass equal to 200 GeV. However, the shape of the 
predicted spectrum, which is the source of the tension with the DAMA data, is effectively 
independent of the DM mass provided it is heavier than about 10 GeV [50j. Therefore the 
curves in Fig. [3] also apply to other DM masses provided we rescale the effective DM-electron 
coupling strength appropriately. 

Our conclusion that electron-DM scattering gives a poor fit to the DAMA spectral data 

■*We also neglect parts of the electron-DM cross section suppressed by the DM velocity. 
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Figure 3: DAMA modulated and unmodulated single hit spectral data and fits of 
electromagnetic elastic DM scattering. The DM particle is assumed to be a fermion of 
mass MjjM = 200 GeV with either scalar or chiral vector {V ± A) couplings to electrons. 
The (red) curve for Fit A is the one that yields the lowest effective value using both 
the modulated and single-hit unmodulated spectral data sets. The (green) curve for Fit B 
corresponds to that with the lowest effective ^^i the modulated dataset alone, excluding 
the lowest energy bin. 



is robust. We find that it continues to hold even if we do not include the lowest modulated 
energy bin and the two lowest unmodulated energy bins in the fit. Furthermore, we have 
also examined other Dirac structures for the couplings between a fermionic DM particle and 
the target electron (relative to the scalar and vector {V ± A) couplings considered above 
and in Ref. [50]), and we find that these do not improve the situation. A scalar DM particle 
scattering elastically off electrons does not appear to work eitherU 



3 General IDM Properties and Nucleon Scattering 

Having found that IDM provides an acceptable fit to the DAMA modulated and unmodulated 
spectral datasets and that other proposals are strongly constrained, we turn now to the 
general properties of potential IDM candidates. Clearly, any such candidate must have an 
inelastic nucleon-scattering cross section that is significantly enhanced relative to its elastic 

^ Inelastic dark matter scattering off electrons might work, although this would likely require an extremely 
large electron scattering cross section. 
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cross section. To meet these criteria, the mass, inelastic sphtting, and inelastic cross section 
must all fall into the appropriate ranges. We now show these requirements can be satisfied 
by IDM particles in the mass range lOOGeV — 5TeV that interact coherently with nuclei 
primarily through the exchange of a massive gauge boson. 

3.1 Inelastic Interactions from a Massive Gauge Boson 

When the DM-nucleon scattering is mediated by a massive gauge boson, the dominance of 
inelastic interactions over elastic arises in a natural way. We consider two types of IDM in 
this context: a Dirac fermion that is split into a pair of nearly degenerate Majorana states 
by a small Majorana mass, and a complex scalar that is split into two real scalars by a small 
holomorphic mass0 

For the fermion case, the models we consider all reduce to 

if D ^tYid^. + tgQ Z'^)^ - M^tfj - Kul {^'^PliI^ + h.c) - ]^mn {^-PriP + h.c.) , (3.1) 

where ■?/' is a Dirac fermion, M ^ ttil^fi, and Z'^ is a massive gauge boson. The Majorana 
masses rriL^R in Eq. (13. ip split the Dirac state ip into a pair of Majorana states \l/i,2- In terms 
of these mass eigenstates, the Lagrangian becomes 

^ D h,iYd^^,-^{M -m+)^,^, (3.2) 

1 - 1 
+ -^f2«7^9^^2 - -(M + m+)^2'^2 

+igQZ'^^>2i^^i 

1 777 _ _ / in^ \ 

+ 2 ^ ^ IF (^27 - ^lYl'^i) + o i^—j , 

where m± = {mi ± m/j)/2. 

From this we see that the dominant gauge boson interaction is strictly off-diagonal, and 
that the mass splitting between the eigenstates is 

5 = M2 — Ml = 2 m+ = niL + rriR. (3.3) 

There is also a residual diagonal coupling of the fermions to the gauge boson, but it is 
suppressed by a power of m_/M -C 1. 

The basic story for the scalar case is very similar. Consider the interactions 

^ D \{d, + igQ Z;)0|2 - M^|0|2 - im2(02 + h.c), (3.4) 

^Another interesting possibility would be spin 1 inelastic dark matter, but models of this type are more 
complicated. 
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where, again, Z'^ is a massive vector boson and we assume 3> rn? with real and 
positive. The holomorphic mass term sphts the real and imaginary components of the 
complex scalar = (0/? + i(j)i)/\/2. In terms of these fields, the Lagrangian becomes 

^ D \{d^nf - \{M' + m^)0?, + \{dh? - \{M' - m')^] (3.5) 

-gQZ'^{ct>id,ct>R - + \g^Q^Z'^Z'^{ct>l + 0?)- 

From this, we see that the single gauge boson interaction with two scalars is strictly off- 
diagonal, coupling (pR exclusively to (pi. The splitting between these mass eigenstates is 

2 

ryt ^ 

5 = VM2 + m2 - VM2 - m2 ~ — . (3.6) 

Let us also emphasize that in both the fermion and scalar IDM cases, it is technically 
natural to have the "Dirac" mass M much larger than the "Majorana" mass m. In the limit 
m ^ 0, both theories have a global IJ{\)y)m symmetry analogous to baryon number in the 
SM, implying that all quantum corrections to m (or m^) are proportional to itself. Indeed, 
within the MSSM the B and L global symmetries keep the real and imaginary components of 
the squarks and sleptons degenerate, while the VEVs of the Higgs complex scalars split their 
components. In Section HI we will construct models that generate such inelastic sphttings in 
simple and natural ways. 



3.2 Nucleon Scattering Rates 

To compute the effective nucleon scattering rates o^^^ relevant for Eq. (12. 2p mediated by a 
massive gauge boson, we concentrate exclusively on vector-vector (VV) interactions. Such 
VV interactions give rise to coherent spin-independent DM scattering off target nuclei [27] . 
Axial-axial (AA) interactions, on the other hand, produce an incoherent spin- dependent 
coupling to target nuclei. Mixed VA and AV interactions can also be neglected because 
they produce effective scattering cross sections suppressed by at least two powers of the DM 
velocity, which is on the order of i; ~ 10~^ in our galactic halo. 

For fermionic IDM arising from couplings of the form given in Eq. (13. ip . the effective 
nucleon cross section cr° needed to compute the interaction rate Eq. (12. 2p is identical to the 
cross section for a Dirac fermion to scatter off the nucleon. Up to small corrections, the 
effect of the inelasticity is completely accounted for by setting the lower velocity cutoff Vmin 
in Eq. (12. ip to the expression given in Eq. (II. ip . Starting with the Lagrangian for vector 
couplings of the SM quarks q and a Dirac fermion to a massive Z' gauge boson (note we 
are using Z' here for generality - the next subsection will restrict attention to the Standard 
Model ZO), 

^^-gg'v Q^'q -gg'tz'^ ^ri^. (3.7) 

the relevant effective nucleon-scattering cross section is [27] 



14 



where Hp^n is the reduced mass of the ■?/'-nucleon system, and 

g^y = 2g^ + g'y, g^=gl + 2g'y. (3.9) 

The couphngs fp and /„ appearing in Eq. ( 12. 2p coincide with gy and gy in the present case. 

For complex scalar IDM, the effective nucleon scattering cross section required to evaluate 
the event rate in Eq. ( 12. ip is identical to the cross section for a single Dirac fermion of the 
same mass as the scalar to scatter off a nucleon. Again the effects of the inelasticity are 
accounted for by modifying the lower velocity cutoff Vmin in Eq. ( 12. ip . With the coupling of 
a complex scalar to a massive Z' given by 

^^-iggtZ'^{<P*d,<P-<pd,<i)*), (3.10) 

the effective nucleon scattering cross section is 

As for the fermionic case, we can identify gy"' with fp and /„ appearing in Eq. (12. 2p . 
3.3 SU{2)l Mediation 

Perhaps the simplest possibility to mediate IDM scattering off nucleons is the gauge 
boson of the SM. Note that the photon is not an option because we assume that the DM 
candidate is neutral and is not an option because the radiatively induced mass splitting 
between charged and neutral components of a multiplet scales as awMw ~ 100 MeV, which 
is much too large for IDM. A neutral particle that couples to the necessarily carries 
hypercharge, so the simplest possibility is that dark matter is a doublet of SU{2)l with 
hypercharge Y = 1/2. 

For a Dirac fermion or a complex scalar SU{2)l doublet with hypercharge 1/2, the 
effective couplings for the neutral components are gy'^^ = 1/2, glr = 1/4 - sin^ew, and 
gy = —1/4. Since siii^ 9w — 0.24, the neutron coupling is much larger than the proton 
coupling. The corresponding cross section is 

^0 ^ 91^1 ^ 7.44 X 10-39 cm2. (3.12) 

271 

To obtain this number, we assumed that the DM mass is much larger than that of the 
neutron so that /i„ ~ m„ = 0.9396 GeV. 

In Figs, m and [5] we show fits to the DAMA modulated data along with constraints from 
CDMS II, CRESST-II, and ZEPLIN-III, for nucleon scattering mediated by Z'^-exchange and 
DM masses of 1080 GeV and 525 GeV, respectively. These particular masses were chosen 
because they are the values that lead to the correct thermal relic density for a Dirac fermion 
(1080 GeV) and a complex scalar doublet (525 GeV) [Hjlll 

^ Using micrOMEGAs v2.2 [55], we find slightly smaller central values for the preferred masses of 
1080 GeV and 525 GeV as compared to Ref. [48], who obtain 1100 GeV and 540 GeV. However, these 
differences are within the margin of error and have a very small effect on the allowed region where the relic 
density scales roughly as floAih^ cc mj-fj^.j. 
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It is intriguing that the inelastic cross section mediated by exchange is very similar to 
the values preferred by our fit to DAMA. This coincidence of scales was observed in Ref. [13] 
in fitting the case of mixed sneutrino DM to the DAMA/Nal dataset. Here we observe 
that this persists even with more detailed fits to the energy spectrum and for more general 
doublet candidates. 

The Z°-mediated nucleon cross section is in fact a little bit too big assuming a local DM 
density of pdm = 0.3 GeV/cm^. However, there is a significant uncertainty in the local DM 
density and lowering its value to poM = 0.15 GeV/cm^, which is also within the allowed 
range, leads to good agreement at the 90% confidence level for fermion doublet DM, and to 
marginal agreement for scalar doublet DM. Another possibility is that the dark matter has 
multiple components [561 [STJ [581 [59j , with the local density of the doublet component giving 
rise to the DAMA signal well below 0.3 GeV/cm^. The fermion or scalar doublet DM could 
also be heavier than 1080 GeV or 525 GeV and its density diluted by a late-time production 
of entropy, allowing for a larger cross section. 

Let us also emphasize that the mass splitting terms in Eqs. (13. ip and (13.41) necessarily 
break f/(l)y. However, such a mass term can be generated after electroweak symmetry 
breaking through operators involving the Higgs field [60]. The exact operator one needs 
then depends on the quantum numbers of the dark matter particle. For S't/(2)2,-doublet 
dark matter with hypercharge y = 1/2, we can introduce the gauge invariant operator 

^ D -^{(j)(f)hh + h.c.) (3.13) 

for a scalar, or 

^ -^{^"-^hh^h.c:) (3.14) 

for a fermion. In order to obtain a splitting 6 ~ 100 keV, we need A ~ {5 /v){m^/v) ~ 10~^, 
or A ~ (v'^/S) ~ 10^ GeV. While without further model structure (such as we will soon 
consider) these are somewhat awkward numbers, we emphasize that the values are technically 
natural. 

Of course one can consider representations aside from an SU{2)l doublet. The next 
simplest possibility is to introduce a complex triplet of SU{2)l with hypercharge Y = 1. In 
order to split the states that couple to the Z'', the triplet would require a higher-dimension 
operator {T^^hihjY in order to be gauge invariant. For scalar DM, this operator should 
be suppressed by a scale A^ ~ {v^/rriTS) ~ (lO^GeV)^, and for fermion DM it should be 
suppressed by a scale A^ ~ ~ (3 x 10^ GeV)^. 

One can obviously keep considering larger representations of SU{2)l, which will cause 
the scale suppressing the smallest gauge invariant splitting operator to decrease even further. 
In extensions to the Standard Model that solve the hierarchy problem, one might naturally 
expect to have operators suppressed by the TeV scale. Large enough representations of 
SU{2)l {e.g. a 5 or a 7) may then naturally have the correct splitting. For a representation 
of dimension N, the DM mass needed to reproduce the right relic density, and hence the 
cross section needed for DAMA, increases roughly as N^'''^. On the other hand, the effective 
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scattering cross section cx^ depends on the choice of hyper charge, and can scale between 1 and 
A^^. This means that thermal dark matter composed of a higher- dimensional representation 
can also agree with direct detection constraints, though the agreement can be better or worse 
according to the direct detection cross sections. 

While it would be interesting to study the fits of larger representations of SU (2) ^ in 
more detail, we find it intriguing that the simplest possibility of an SU{2)l doublet works 
reasonably well for explaining the DAMA data. If one is to take this model of dark matter 
seriously, then, the main question is whether the numbers and scales cited above in Eqs. (13.131) 
and fl3.14p . though technically natural, have a reasonable origin in models. That is, suppose 
DAMA has indeed discovered dark matter. What would be a reasonable interpretation of 
this result? In Section HI we will consider several possibilities for explaining the physical 
origin of these operators and their coefficients. 
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Figure 4: DAMA allowed region and constraints for scattering through Z^'-exchange for a 
fermion SU{2)l doublet with a mass of 1080 GeV, a value that yields the correct thermal 
relic density. We consider values of the DM escape velocity of Vesc = 500 km/s and 600 km/s 
and local DM densities 0.15 GeV/cm^ and 0.3 GeV/cm^. 
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Mdm = 525 GcV, Vese = 500 km/s, p = 0.3 GeV/cm' Mdm = 525 GeV, = 600 km/s, p = 0.3 GeV/cm' 
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Figure 5: DAMA allowed region and constraints for scattering through Z^-exchange for a 
scalar SU{2)i doublet with a mass of 525 GeV, a value that yields the correct thermal relic 
density. We consider values of the DM escape velocity of Vesc = 500 km/s and 600 km/ s and 
local DM densities 0.15 GeV/cm^ and 0.3 GeV/cm^. 
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3.4 U{1)^ Mediation 



Inelastic DM scattering can also be mediated by a massive gauge boson not existing in the 
Standard Model. Two distinct possibilities that can yield the right cross sections are a heavy 
(TeV-scale) gauge boson with order unity couplings to Standard Model fields, and a light 
{Mz' -C Mzo) hidden gauge boson with highly suppressed couplings to the SM. If the hidden 
gauge symmetry is an Abelian U{l)x, small couplings arise in a natural way from kinetic 
mixing with hypercharge. We consider here both the heavy and light exotic gauge boson 
cases, focusing on an Abelian U{l)x gauge symmetry for simplicity. 

3.4.1 A Heavy Visible f/(l)^ 

A heavy U{l)x Z' gauge boson with order unity couplings to the SM can confiict with 
phenomenological bounds on the mass of this new state. Collider and other bounds place 
lower limits on the Z' mass, whereas a very heavy gauge boson generate a nuclear scattering 
cross section that is too small to account for DAMA, yielding some tension in this scenario. 

Precision measurements at LEP imply that the bounds on lepton couplings for a given 
Z' mass are generally stronger than those for quarks. For example, a {B — L) gauge 
boson would not satisfy phenomenological bounds and allow for a DAMA signal without 
unreasonably large couplings to the DM particle or an extremely small gauge coupling. 
Satisfying phenomenological constraints, even for a gauge coupling to the SM as small as 

~ 0.4, requires a gauge boson mass Mz' ^ 2.5 TeV. To achieve the DAMA signal would 
then require the effective coupling to dark matter to be greater than QxXdm ^ 10, where 
xdm is the dark matter charge. Lighter (B — L) gauge bosons with perturbative couplings 
to the DM {gx xdm ^ 1) are possible at the expense of making the gauge coupling gx < 0.05 
while keeping x^m ^ 20. Such a hierarchy of charges seems contrived. 

Models with smaller coupings to leptons relative to quarks and the DM are more rea- 
sonable. In Fig. [6] we show the values of the U{l)x charges of the leptons {xe and xl for 

and L) and the DM particle (xdm) for allowed points from a scan over heavy U{l)x 
models. These points satisfy both the phenomenological constraints on a heavy Abelian 
gauge boson and generate a reasonably large nuclear scattering cross section. We assume 
fiavor-independent gauge charges (consistent with Yukawa couplings) for simplicity and to 
minimize fiavor- violation, and take a two-Higgs doublet model for generality. Anomaly 
cancellation can be satisfied by adding exotic fermions. We have focused on charges such 
that gxXi < 1 to ensure weak coupling. 

The phenomenological bounds we apply in generating Fig. [6] are the direct search bounds 
from the Tevatron in the di-lepton [61j and di-top [62] channels and the limits on contact 
interactions from LEP II [63]. In applying the Tevatron bounds, we assume the Z' decays 
entirely into SM final states. We do not explicitly compute the effects of the new U{l)x on 
precision electroweak observables, but demand that the Z^ — Z' mixing angle be less than 
(^mix < 3 X 10"^ [Ml ESI EU EZl [HH [69], which approximately captures the constraints from 
these observables [66]. In computing the mixing angle, we assume that tan/3 = 10. For the 
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Figure 6: Values of the lepton and dark matter charges under U{l)x consistent with the 
phenomenological constraints hsted in the text that also generate an adequately large value 
of the effective nucleon scattering cross section o'e//; as a function of the Z' mass Mz'. The 
charges of the quark and Higgs fields are also scanned over. 



nucleon scattering cross section, we demand that the quantity 



be larger than 3 x 10 cm? or 3 x 10 cm?. This definition and these lower limit values 
are motivated by our fits to IDM in Section [2|, where we assumed fp = Qy = fn = Qv 
Eq. (D- 

We conclude from Fig. [6] that a Z' with a mass in the range of 400 GeV to several TeV can 
generate an adequately large nucleon scattering cross section to account for DAMA while not 
running afoul of the existing bounds. Such a Z' must generally be somewhat leptophobic, 
and have significantly large couplings to quarks and the DM state. However, we note these 
couplings are so large that the corresponding gauge coupling generally encounters a Landau 
pole well below the GUT scale, even without including the effects of possible exotics required 
for anomaly cancelation (which would only make this problem worse). Thus, many popular 
U{l)x models such as {B — L) and the exotic f/(l)'s motivated by Eq [70j do not work if 
a unification relation is imposed on the exotic gauge coupling. Furthermore, although this 
heavy Z' scenario is a possibility, the natural connection to the thermal relic abundance that 
was present with an electroweak cross section is in general lost. 



3.4.2 A Light Hidden f/(l)^ 

A second possibility for an exotic Abelian gauge symmetry mediating IDM scattering at 
DAMA consists of a relatively light hidden U{l)x that couples only weakly to the SM. Very 
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small SM couplings arise naturally if the hidden sector couples to the SM through kinetic 
mixing of the U{l)x with hypercharge. A coupling [7l] 

if (3.16) 

where B^^i, and X^j, are the otherwise canonically normalized U{1)y and U{l)x field strengths, 
can arise from loops of heavy states that are charged under both gauge groups. This leads 
to typical values of e ~ 10~^ — 10~^ [72]. Making a field redefinition to eliminate the kinetic 
mixing term, the SM matter fields acquire effective charges under the U{l)x- If the DM 
states couple directly to the U{l)x but not the SM gauge groups, the exotic gauge boson can 
mediate DM scattering off nuclei. 

We are interested in the case where the exotic Z' gauge boson is light, Mz' <^ Mz- In 
this limit the U{l)x can be treated as mixing with U{l)em, and the induced charges are 

^ D e coseweQZ'jYf + 0{eMl,/Ml), (3.17) 

where / represents a SM fermion with electric charge Q. Since the light Z' couples to electric 
charge in the visible sector, it will mediate scattering of the DM only with protons. Applying 
Eq. 03.81) or Eq. fl3.1ip . the corresponding cross section is 

where QxXdm is the charge of the DM particle under the U{l)x- Thus, for gauge boson 
masses of a few GeV and values of e towards the lower end of the reasonable range, the 
DAMA signal can arise from IDM scattering off nuclei mediated by a hidden Z' . 

Phenomenological bounds on a light hidden U{l)x have been studied in Refs. [66 | [73 | [7 ^ 
[75l [76] . The strongest constraints on a hidden U{l)x gauge boson with mass of a few GeV, 
coupling to the SM through kinetic mixing with electromagnetism, comes from measurements 
of the magnetic dipole moments of the electron and muon. For masses larger than about a 
GeV, values of the kinetic mixing parameter e less than 10~^ are acceptable . 



If the mass of the dark matter particle is much larger than the mass of the light U{l)x 
gauge boson, a potential fermionic IDM candidate and its slightly heavier partner will both 
annhilate very efficiently into gauge boson pairs. The annihilation cross section for this 
process is [77] 

Note that this cross section scales as ~ g'^/M'jji^j, just as for electroweak gauge boson 
exchange. For reference, the necessary thermal DM relic density is approximately ^dm — 
2 (3 X 10"^'^ cm •^/s)/ where the additional factor of two accounts for the fact that 

two different states are annihilating [TSl [79]. In addition, there are contributions to the 
annihilation cross sections from the s-channel exchange of U{l)x gauge bosons decaying into 
pairs of light matter fields in the U{l)x sector. 
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These contributions depend on the precise field content of the theory, but they will 
typically be on the same order as for annihilations into light gauge boson pairs. In general, 
the correct thermal relic density of a fermionic IDM state in this scenario will be obtained 
for a dark matter mass on the order of several hundreds of GeV. For a complex scalar IDM 
state, there is a moderate additional velocity suppression of the annihilation rate, and a 
somewhat lighter mass (but still on the order of a few hundred GeV) will yield the correct 
thermal relic density [73] . 

It follows from the discussion above and Eq. flS.lSp that the we can choose parameters to 
give the correct thermal abundance while simultaneously accounting for the DAMA signal. 
Recall, however, that in the case of electroweak interactions this came out automatically for 
the known gauge boson masses and gauge coupling. 

4 Models of Inelastic Dark Matter 

Our guiding assumption in investigating candidates for IDM is that the interaction between 
the DM state and the target nucleus is mediated by a massive gauge boson. Though it 
is technically natural for the mass splitting operator coefficients to be small because they 
violate a global U{1)dm symmetry, the primary challenge is to generate a small inelastic mass 
splitting in a reasonable way. With this in mind, we introduce a f/(l)DAf-breaking spurion 0. 
Depending on the relative charge of the dark matter field and 0, the mass splitting operator 
will be naturally suppressed by a factor (^f^)". This operator could be used to generate a 
~ 100 keV splitting for a judicious choice of n and < (p >. 

This is not terribly satisfying, however, as it simply parameterizes our ignorance about 
how the symmetry is broken. We would also like to understand the physics underlying the 
inelastic splitting and, given that the masses generally hover around the electroweak scale, 
fit it together with possible solutions to the hierarchy problem. The goal of this section is 
to come up with simple models that can give rise to a splitting of the correct size in the 
context of a solution to the hierarchy problem without any large tuning of parameters. Of 
course, in the end aesthetic criteria are subjective, so we view this section as a compilation 
of interesting ways to generate splittings of the right size for fermions or scalars. 

4.1 Models Mediated by the SM 

As we saw in Section |3l dark matter charged under SU{2)l seems particularly promising 
because the Z°-exchange cross section is roughly the correct size to account for the DAMA 
signal when one chooses a dark matter mass that gives the right thermal relic abundance. 
This is especially interesting for SU{2)l doublet dark matter, which is the simplest possible 
representation. However, the size of the mass splitting is then somewhat of a mystery, with 
scalar dark matter requiring the coefficient of the splitting operator Eq. (13.131) to be A ~ 10~^, 
and fermion dark matter requiring the splitting operator Eq. (I3.14p to be suppressed by the 
mass scale A ~ 10^ GeV. Both splittings are smaller than what one would naively expect 



23 



for an effective field theory valid below the TeV scale if U{1)dm were strongly broken. 

One simple possibility is that the effects of U{1)dm breaking are sequestered in some 
way. This could happen, for example, if the splitting operators are generated by integrating 
out a singlet S through an operator DS*h, where h is the SM Higgs field and D is the dark 
matter doublet. In this case, U{1)dm breaking could be communicated through the singlet. 
If the singlet is very heavy, or only couples weakly to the doublet D, the mass splitting is 
suppressed. Note that if there are multiple trilinear couplings DS*h and DSh, we can simply 
call the linear combination that couples to the doublet S* and define the U{1)£)m symmetry 
so as to respect this trilinear coupling. 

There are then two ways that f/(l)z)Af breaking could be communicated. The first 
possibility is that most of the singlet mass is U{1)dm preserving, with a small U{1)£)m 
breaking piece. For a scalar with the potential 

^ D -ml\S\^ - (^^S^ + fDS*h + h.c^ , (4.1) 

integrating out the singlet generates the operator ^-^DDhh. For a fermion with Lagrangian 

^ D -msSS - (^^S^S + XDSh + /i.e.) , (4.2) 

integrating out the singlet generates the operator ^J^^ D'^Dhh. If there is a natural hierarchy 
between ms and mg, or if the couplings A or / are naturally small, one could obtain a splitting 
of the correct size. 

The second possibility is that the singlet has a very large U{1)dm breaking mass rris. 
This would, for example, be the only option if the singlet is a real scalar or a Majorana 
fermion. Integrating out the singlet would then generate the scalar operator -^DDhh or the 

fermion operator ^D'^Dhh. This is completely analogous to the way small neutrino masses 
are induced in the conventional seesaw mechanismj^ A splitting of the right size could again 
be obtained for small couplings or if 7715 is naturally identified with an intermediate scale. 

Only a handful of concrete models exist for IDM whose nucleon scattering is mediated by 
the Z^. Ref. [TTj proposed a model of left-handed sneutrino dark matter in which the inelastic 
mass splitting is generated through mixing with a scalar singlet right-handed sneutrino. 
The f/(l)_DM violating (and lepton-number violating) mass arises through a SUSY breaking 
operator j^X^X^XN^N, and the size of the splitting is naturally related to an intermediate 

scale. For related models, see Refs. [501 EI]- A second possibility for IDM are the Dirac 
neutralinos that arise in f/(l)/{-symmetric SUSY scenarios [121 [IT]. The small Majorana 
mass splitting would then be related to a small amount of f/(l)/j breaking. 

We present below several other models for SU{2)l doublet IDM that make use of a 
singlet to communicate U{1)dm breaking. Two of these models are based on a warped extra 



®In fact, it may be possible that the same right-handed neutrino scale enters both the neutrino mass and 
the dark matter splitting operators - this would require that the lepton doublets come with an additional 
suppression factor, but this may be natural in a model of flavor physics. 
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dimension [82], and illustrate some of the ways an inelastic splitting can emerge in this 
context. We also present a supersymmetric model. 



4.1.1 Warped Fermion Model 

We begin with a model of fermion SU{2)l doublet dark matter and attempt to explain 
how the scale of splitting 6 ~ 100 keV can emerge without any large hierarchy of input 
parameters. This requires an explanation for the scale suppressing the splitting operator 
Eq. fl3.14p . A ~ 10^ GeV, which could represent the mass of a singlet field that has been 
integrated out. One way that this intermediate scale mass could emerge naturally is if it is 
equal to the Planck scale times an exponential suppression factor. As we will see below, it is 
straightforward to realize this possibility in the context of a 5D Randall- Sundrum model j82j . 
giving the added bonus of combining a natural dark matter model with a solution to the 
hierarchy problem. Our model is similar to the models of Refs. |83l [H], which realize the 
seesaw mechanism in warped geometry. 

In particular, we consider AdS^ compactified on 5*1 /Z2 with metric 

ds^ = e-^^^y^r^^^dx^'dx" - dy^, (4.3) 

where —nR < y < ttR. If the dark matter derives from a vector-like SU (2) l fermion doublet 
D = {DljDji)^ localized on the TeV brane, its mass is naturally of order the TeV scale. 
However, since we also expect the Higgs doublet to be localized to the TeV brane, we need 
to forbid the TeV brane localized operator -j^ — D'^Dhh which would generate too large of a 
splitting. We therefore impose a U{1)dm symmetry under which D is charged to forbid the 
splitting operator, and assume that this symmetry is broken only on the UV brane. 

In order to communicate the breaking of U{1)dm, we introduce a bulk fermion singlet 
S = {Sl,Sr)^ with (+, +) boundary conditions for Sl- We include a ?7(l)DM-breaking 
Majorana mass on the UV brane, along with a [/(l)i:)M-pi'eserving bulk mass. The singlet 
action is taken to be 



S = d'^x / dyy/^ 



iS-i^^DMS + cke{y)SS - 6{y) ( ^SISl + h.c 



S{y - 7tR) (^XDnSLh + h.c^j 



(4.4) 



where only the left-handed component of the singlet can have brane couplings because of 
the choice of boundary conditions. Here, S'^ = C'y^S* where C is the 5D charge conjugation 
operator. The sign of the bulk mass parameter c has been chosen to agree with the 
convention that c > 1/2 localizes a zero mode towards the UV brane. Note that while the 
boundary singlet mass explicitly violates the f/(l)£)M symmetry, there remains an unbroken 
Z2 subgroup under which D and S are odd ensuring that the lightest of these fermions is 
stable. 

To see the effect of the boundary terms on the fermion masses, it is easiest to first expand 
the singlet in a basis that diagonalizes the KK modes without the boundary terms. In this 
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basis, the communication of U{1)dm breaking is dominated by the chiral zero mode which 
picks up a large Majorana mass. Since the remaining modes acquire Dirac KK mass terms, 
we can truncate the KK tower while still capturing the dominant contribution. In particular, 
one can expand 



where the wavefunctions /l „ solve the bulk equations of motion 

{dy±ck)fL Jy) = Tm„e^l^l/fl_„(?/). 



(4.5) 



(4.6) 



Imposing (+, +) boundary conditions for Sl, the solutions are 



My\/2 



,m 



(4.7) 



where the masses m„ can be determined from the condition that /i?,n(0) = 0, and the 
normalization factors are obtained from 



1 



In particular, there is a massless chiral zero mode with 



(4.8) 



I _ g-(2c-l)7rfc_R, 



(4.9) 



The other KK masses and normalization factors can be approximated for m.„ « k and 
fci? » 1 as [85J 



TTlr. 



(4.10) 



assuming c < 1/2, and 
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(4.11) 



After performing this decomposition, the 4D fermion mass matrix is not yet diagonal due 
to the brane localized mass terms. The singlet mass matrix is 
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(4.12) 
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where the Majorana masses Amn are given in terms of the wavefunctions on the UV brane 

A^n ^ ^/^ „^(0)/^,„(0). (4.13) 

In particular, the zero mode picks up a mass 

_ duvkjc- 1/2) 

The couphngs to the canonically normahzed doublet and Higgs field on the TeV brane 

are 

oo 

^ D - J2CnDRS2h + h.c. (4.15) 

n=0 

where 

C„ = -^f^^,(nR) (4.16) 

is determined from the wavefunction overlap on the TeV brane. Assuming that we choose c 
such that Aqq » nin, we can to a good approximation simply integrate out the heavy chiral 
mode. This generates the splitting operator 

S£ D -^D%DRhh + h.c. (4.17) 



2^00 
2d 



e-^(-~m->'RD'j^DRhh + h.c. (4.18) 

uv 

Choosing the natural values ~ -^J-^ and duv ~ 2, we need c ~ 0.13 to obtain a mass 
splitting on the order of 5 ~ 100 keV. It is straightforward to check that including the KK 
modes gives a negligible (~ 10^^) correction to these estimates (see Appendix . This is 
due to the fact that, unlike the former zero mode, the states in the KK tower are mostly 
Dirac and are very poor at communicating U{1)om breaking. 

Arranging the model parameters to produce a mass splitting on the order of 5 ~ 100 keV, 
our model has all the ingredients needed to account for the DAMA annual modulation result. 
The lightest stable fermion state is an almost pure SU{2) l doublet Majorana fermion with 
a slightly heaver inelastic partner, and is stable on account of the unbroken Z2 subgroup of 
U{1)dm- If the corresponding fermion mass is ~ 1.1 TeV, it will the yield the correct relic 
density from thermal freeze-out, in addition to the correct nucleon scattering cross section 
and mass splitting as was shown in Section [31 



4.1.2 Warped Scalar Model 

We now turn to a model of scalar SU{2)l doublet dark matter. To sequester U{1)dm 
breaking we again consider a warped 5D setup with the metric of Eq. (14. 3p . and assume that 
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a complex scalar doublet D = -^{011 + iDj) is localized to the TeV brane in addition to the 
Higgs field. As in the previous model, we need to forbid the TeV brane-localized operator 
DDhh, which would split the masses of Dr and Dj by ~ with an 0(1) coupling. We 
again assume a U{1)om symmetry which is broken only on the UV brane, while preserving 
a Z2 subgroup to ensure the stability of the doublet DM. 

In order to communicate this U{1)£,m symmetry breaking to the DM doublet, we proceed 
as in the fermion model and introduce a bulk (scalar) singlet. Note that we cannot allow 
the doublet to directly couple to the UV symmetry breaking, even though a wavefunction 
localized to the TeV brane can be sufficiently suppressed in the UV, because the UV mass 
would then break hypercharge (since a neutral dark matter doublet must carry hypercharge) 
at too high a scale. Instead we assume the IR-brane-localized doublet mixes with a bulk 
singlet, as in the fermion model of the previous section. 

The singlet can then couple in symmetry-breaking operators on the UV brane and 
communicate this breaking to the doublet on the TeV brane. However, we note a critical 
difference between 5D scalars and fermions - in general bulk scalars do not possess a zero 
mode. This means that the only option is to communicate U{1)dai breaking through singlet 
KK modes, which naturally peak away from the UV and tend to pick up small U{l)DAr 
breaking masses (in addition to their [/(l)£)M-preserving masses). Therefore this model 
does not suppress mass splitting via a seesaw mechanism. Instead, the real and imaginary 
components of the doublet will be split by their mixing to the split components of the bulk 
scalar KK modes. 

In fact, bulk KK modes suppress the communication of U{1)£,m symmetry breaking to 
the doublet extremely effectively - so much so that the effective UV scale cannot be much 
higher than ~ 100 — 1000 TeV or the U{1)dm breaking would effectively decouple and give 
too small a mass splitting to the dark matter candidate. In the context of an RS solution 
to the hierarchy problem, this requires a setup with either a third brane|^ or an additional 
warped dimension such as was considered in Ref. [88]. Alternatively, one can simply treat 
the smaller warp factor as giving a solution to the flavor hierarchy problem, such as in the 
little RS scenario of Refs. [89l [90] . The model of this section can be thought of as a two-brane 
effective theory describing any of these situations, generated by integrating out physics above 
the effective UV scale. 

We thus start by introducing a complex bulk singlet S with (+, +) (Neumann) boundary 
conditions. In general, S has a bulk mass and two [/(l)DM-preserving brane mass terms, 
as well as a f/(l)£)M-breaking mass term on the UV brane. To simplify the calculation we 
assume that the [/(l)i3M-preserving brane mass terms are negligible. We consider the singlet 
action 



S = j d^x j dy^ ^dMS^d'^'S - m^S*S - 6{y) + h.a 

- 6{y - ttR) (a e^"'^^ DS*h + h.c.)] , (4.19) 



^Detailed considerations of the three-brane setup, such as ensuring its stability, are beyond the scope of 
this work. Related references can be found in, e.g., Refs. [86l [87], 
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where m = -Jak is the f/(l)DAf -preserving bulk mass of muv is the [/(1)£)m- violating 
UV brane mass, and A is the coupling to the dark matter field D. The exponential factor is 
because D and h are taken to be canonically normalized 4D fields, so A ~ ^/k has units of 
Y^mass. 

The KK decomposition of S is 



S{x^, y) = -j= J2 S-{x^)Uy), (4.20) 



n 



where the functions /„(y) that solve the bulk equations of motion are 



/n(l/) = ^ Jo.(^e'\'\) + h^{mn)Y^{^e'\y\) . (4.21) 



2k\y\ 



Here, a = a/4 + a, and the normalization factor Nn can be approximated in the limit m„ << 
k and kR >> 1 as 



Nn-x^ (4-22) 

The masses m„ and the functions ba are determined from the boundary conditions. In the 
f^(l) DA/- symmetric limit {muv = 0), the boundary conditions of Sr and Sj are identical, so 
there is a pair of degenerate states at each KK level. When U{1)dm is broken {muv 7^ 0), 
Sr and Sj have different boundary conditions and the resulting KK masses are split. 

In the presence of the f7(l)DM-violating mass, variation of the boundary action deter- 
mines the boundary conditions to be 

dySR^muvSn = 0\y=o (4.23) 

It is straightforward to solve numerically for the KK mode mass splitting in the presence of 
these boundary conditions. This mass splitting in the singlet sector is then communicated 
to the DM doublet through the couplings on the TeV brane. 

It is easier however to proceed as in the last section and first expand the KK states in the 
basis without the splitting, treating the f/(l)£)Af-violating mass as a perturbation. Imposing 
the simpler boundary conditions dyS = 0|j,=o,7ri? leads to the two conditions 

9 T ( _|_ UhL J' (mn \ 

b (m ) = - ^ ' k -Jgy k ) u 24) 

^ia\ I, }^ k ^aV k f 



" k '' k k l_ 

oVfc I k k ' 



(4.25) 



which yield the approximate KK spectrum in the limits « k and kR >> 1 [85j, 

m„ ^ + I - ^ j TT A; e'''^^. (4.26) 
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In this basis, the masses in the singlet sector can be written as 

^12 
^2 I A 

^22 



^ D -1(3}, SI ... ) I ml + Al 
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^21 '"'2 
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~=^{s] s] ... ) 
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^21 



12 

ml - Aj^ 




(4.27) 



where 



(4.28) 



The mass sphtting in the singlet sector at each KK level is then determined by the 
difference in the eigenvalues of these two mass matrices. As long as we choose parameters 

the mixing is small and the eigenvalues are simply ± A^ 



such that Amn << rnin-, 

O(^) corrections. Expanding the wavefunction profiles in Eq. fl4.2ip for m„ << k at y 
the mass splitting at the n**^ KK level is then approximately 
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(4.29) 



The mass splitting in the singlet sector is transmitted to the DM doublet through their 
coupling on the TeV brane. After the Higgs acquires a VEV, the mass matrix is 
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(4.30) 



where 



Cn 



-2-KkR 



■UttR). 



(4.31) 



It is straightforward to diagonalize this mass matrix at a given level of truncation of the 
KK tower. For example, including just the first KK state leads to a mass squared splitting 
to leading order in An of 



Ami, 



A^ 
^11 



mT 



m 



D\ 



(4.32) 
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Note that this goes to zero if either the UV brane mass or TeV brane couphng turns off, 
just as one would expect. We can further expand this for small v « rrin, and write the 
contribution from the n**^ KK level more generally as 



AniD ~ Am„^^ 1 + , (4.33) 

where Am„ is the singlet mass splitting at the nth KK level, given approximately in Eq. fl4.29p . 
Choosing the values k ~ 500 TeV, R ~ 2.1/ k, rriuy ~ 2/c, a ~ 0.1, A ~ 2.5v^, and 
m/) ~ 525 GeV, for example, gives a contribution to the DM mass splitting from the first 
KK mode of ~ 5 keV0 Note that a splitting in the singlet sector that is not too small 
requires the warp factor kR not be too large. 

However, when trying to sum these KK contributions, the sum apparently diverges. This 
is counter- intuitive, since we would expect heavy modes to decouple. Here we might expect 
this decoupling to happen since the contribution to the splitting appears to be smaller for 
heavier KK modes due to the ~ 1/n^ suppression in Eq. ( 14.33^ . However, this is compensated 
for by the fact that the singlet splitting Am^ is increasing for higher KK modes due to a 
larger UV brane overlap, and scales roughly as ~ n^" > n^. 

The issue then is the behavior of C„ for higher KK modes. In the naive limit of an 
infinitely thin brane, the higher KK modes would appear to couple to the TeV brane with 
roughly equal strength, and the magnitude of Cn would not significantly change as one 
increases n. This behavior was seen, e.g., in Refs. [91] [92] . where the IR brane coupling of 
KK SM gauge fields was studied and found to be universal. This would lead to a contribution 
to the DM mass splitting that increases approximately linearly with n for small values of 
the bulk mass a. 

However, this result is unphysical since it does not take into account the thickness of 
the IR brane, A ~ A"^, where A is the 5D cutoff scale and the thickness is given in y- 
coordinate space. We expect A ~ (10 — 100) A; based on naive dimensional analysis. Note 
that the physical thickness of the brane at y ~ ttR is redshifted to Aphys ~ h'^e^^^. The 
effective 4D coupling C^^^ is then better approximated by integrating the wavefunctions 
over this thickness. While not terribly important for lower KK modes, higher KK modes 
rapidly oscillate over this region and the effective coupling is suppressed due to a cancelation 
between opposite phases. 

While the exact numerics depend on the details of the Higgs and DM doublet wavefunc- 
tion profiles over the brane thickness, we can understand this decoupling by assuming the 
"flat" profiles 

h^D{y) ^ ^e^\ (4.34) 



^A 

and concretely taking A^^ ~ 10 k. The effective coupling C^^^ is then given by 

C:ff^^==- dye-^'^yUy). (4.35) 



^°Note that while kR 2.1 is not much larger than 1 as was assumed in some of the approximations above, 
we find that the analytical formulae still give a reasonable approximation to the full numerical results. 
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With these choices, only the first ~ A/A; ~ 10 KK modes make a significant contribution to 
the DM mass sphtting. Taking the same parameters as before, we obtain a mass sphtting 
of ~ 70 keV when summing over the first 10 KK modes. Note that this decouphng behavior 
itself is quite robust in that it will happen for any reasonable choice of Higgs and DM doublet 
profiles. 

For higher values of n, the integrand in Eq. (14.351) is rapidly oscillating over the brane 
thickness A. This happens when the KK modes have enough 5D momentum such that 
many wavelengths A„ fit inside the physical thickness of the brane. We then expect a phase 
cancelation up to terms suppressed by about ~ \n/ ^phys ~ l/'Ti (see Appendix IB] for a more 
detailed derivation). For small bulk mass a, Am„ in Eq. (I4.29p scales roughly as ~ n^, so 
the mass splitting in Eq. (14.331) now scales as ~ 1/n. Summing the contributions from these 
higher KK modes up to n ~ l^^-^kR (corresponding to modes with momentum around the 
5D cutoff scale) then gives an additional contribution which is enhanced by ~ log(e'^'^^) ~ 6 
relative to the contribution from a single lower mode. Estimating the overall amplitude as 
~ 1/10 of the sum of the first 10 modes, these higher KK modes can then give an 0(1) 
correction to the mass splitting, and one can obtain a splitting of the desired size. 



4.1.3 A Supersymmetric Candidate 

The mechanism described above for generating a splitting, namely mixing the DM doublet 
with a singlet that has a [/(l)£)A'/-breaking mass, can apply in other contexts as well. We 
next consider the application of this mechanism in the context of low-energy supersymmetry. 
In this case, a natural small suppression can arise for example in large tan j3 scenarios when 
mixing the singlet with the doublet dark matter candidate through the VEV of Hd. 

An example of a model that exploits this suppression consists of the MSSM augmented 
by a vector pair of SU{2)l doublet chiral superfields D and D'^ with Y = ±1/2. Hypercharge 
and holomorphy then allow the superpotential operator DD'^, but not DD or D'^D'^, and an 
accidental U{1)dm global symmetry can arise. This symmetry must be broken to generate 
an inelastic splitting. To do so, we gauge a related U{l)z symmetry and break it by Higgsing. 
A simple superpotential that realizes the above symmetries is 

W D XNHu-Hd + X'SHd-D + ^NS^ + CNDD". (4.36) 

Here, and S are SM singlets that carry non-zero charges under the gauged U{l)z group. 

If the A^ field develops a VEV induced by supersymmetry breaking soft terms, A^ 
(N) ~ TeV, a supersymmetric mass for S will be generated^ To get an approximate 
picture of what this does to the masses of the fermions in the model, we can integrate out 
the S superfield in the supersymmetric limit. This yields 

Weff D X {N) Hu-Ha + C {N) DD^ - — — {Ha-Df. (4.37) 



'^-'^The size of the N VEV can be set by SUSY breaking so that it is naturally on the order of a TeV. The 
A'^ scalar will get a mass upon expanding the scalar potential, while the N fermion will develop a mass by 
mixing with the U{l)z gaugino. 
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The last term is the desired mass sphtting operator. It receives a suppression from large 
tan/3, our choice that A' < 1, and from an assumed small hierarchy between (N) and the 
electroweak scale. A more careful analysis of mixing in the SM-neutral fermion sector after 
symmetry breaking shows that the fermion mass eigenstates consist of an almost pure singlet 
and a nearly degenerate pair of Majorana states that derive almost entirely from the doublets. 
The resulting SM-neutral fermion masses are 

M,, =C (iV)+5±, M, = ^ {N) + 6+-6., (4.38) 

2 

with 

,,^±^fl±§lV'. (4.39) 



2M, V Ms, 

Taking tan/3 = 30, M, = 3000 GeV, = 1000 GeV, and A' = 0.1, we obtain a mass 
splitting of 5+ — 5_ ~ 130 keV. 

The neutral fermion components of D and D'^ can therefore yield IDM provided these 
fields are stable. An unbroken Z2 discrete symmetry is the minimal possibility to ensure 
this. If this symmetry is i?-parity, D and D'^ as well as the and S superfields must all be 
even. On the other hand, with a new Z2 it is possible for D, D'^, and S to be odd, with A^ 
and the Higgs fields even. The inclusion of such a new discrete symmetry is well-motivated 
in gauge-mediated models with a light gravitino. 

The form of the superpotential in Eq. fl4.36l) can be enforced by the U{l)z charges 

[Sl = Zs, [N], = -2zs, [Dl = Zd, [D% = 2zs-Zd, (4.40) 
[Hd\z = —Zg — Zd, [Hy\z = 3zs + Zd- 

With these charges, the dangerous operator [Hu-d'^]z = 5zs is forbidden provided Zg 7^ 0, 
as is the bare /i term operator [Hu-Hd]z = 2zs. Note also that if Zd = Zg, a Z2 that 
stabilizes both D and D'^ arises automatically as an unbroken discrete subgroup of the U{l)z 
gauge symmetry. As it stands, this theory has mixed S'f/(2)| x U{l)z, U{1)y x U{1)z, and 
U{1)1 X U{1)y anomalies when Zg 7^ 0, implying that SM-charged exotics must also be 
present in the theory. In fact, such exotics are necessary in any MSSM gauge extension that 
forbids a bare /i-term [HSIIHZI- These exotics need not interfere with the dynamics discussed 
here. 

It is more challenging to obtain an acceptable scalar inelastic dark matter splitting from 
this model. The scalar components of D and D'^ will be stable if these superfields (and 5*) 
are odd under i?-parity or if we impose a new Z2 discrete symmetry. In the latter case, it is 
necessary to tune the soft masses such that the fermion components of D and are heavier 
than the scalars. At the supersymmetric level, a scalar mass splitting arises from the F-term 
potential for Hd derived from the effective superpotential of Eq. (14.371) 



A 



/2 



Hd-DD + n.ffH^ 



(4.41) 



D Md DD + h.c, 
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with Md and Ms as in Eq. (14.381) . and /Xe// = A (A^). Relative to the fermion sphtting of 
Eq. (I4.39p . this operator is suppressed by only a single power of cos/3 ~ 1/ tan/? at large 
tan (3. If generic soft supersymmetry breaking operators are also included, there arises a 
further scalar mass splitting with no cos (3 suppression at all. This contribution can be 
thought of as coming from an operator of the form 

where is the trilinear soft parameter corresponding to the N S"^ superpotential operator. 
Numerically, we find that the operators of Eq. (14.411) and (I4.42p produce too large of an 
inelastic scalar mass splitting unless A' < 0.01 and there is an additional small hierarchy 
between fJ^eff, and M^, Mg. Note that even when the scalar splitting is too large, the 
radiative corrections to the fermion mass splitting are still safely small. 

Scattering off nuclei by the IDM candidates that arise in this model will be mediated 
primarily by the SM . The massive U{l)z gauge boson can also contribute to nuclear 
scattering, but the effect will be suppressed by its larger mass. If the new states in the 
U{l)z sector are somewhat heavier than the IDM candidate, the thermal relic density will 
be determined primarily by electroweak interactions, and our estimates from Section [3] for 
the relic density carry through. An interesting additional possibility arises when the doublet 
IDM state is stabilized by a new Z2 symmetry, rather than i?-parity. In this case, the lightest 
superpartner will provide a second contribution to the dark matter. Such multi-component 
DM scenarios have been considered in a number of recent works [FTf ESI ESj- As long 
as any additional DM component has a small scattering cross section off nuclei, it will not 
ruin the IDM explanation for the DAMA signal provided the IDM component makes up a 
significant fraction of the DM relic abundance. 

4.2 Models Mediated by an Exotic Z' Gauge Boson 

In Section [3] we showed that the DAMA signal can arise from the scattering of IDM off iodine 
nuclei mediated by an exotic U{l)x gauge boson. We distinguish two possibilities for such 
a gauge boson, one in which it is heavier than the SM and couples directly to the SM 
and the other in which it is much lighter than the Z^ and couples only weakly to the visible 
sector, such as through a small kinetic mixing with electromagnetism. We construct here 
models for IDM that realize both possibilities, though we will see the latter involves mass 
scales that are more contrived. 

4.2.1 Heavy f/(l)^ Models 

The Z°-mediated candidates for IDM presented above can all be adapted to models in 
which the scattering off nuclei is mediated by the exchange of a heavy U{l)x gauge boson. 
In each case, the SU{2)l doublet states in our previous models are replaced by a pair of 
states with vector-like charges under a U{l)x gauge symmetry. For the analog of our SUSY 
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model, the U{l)x gauge symmetry should be broken by a pair of fields such that they obtain 
hierarchically different expectation values from the dynamics of the potential. 

In contrast to the Z'^-mediated models presented above, however, the thermal rehc density 
in heavy U (l)^ models is more model dependent. If the DM particle is lighter than the exotic 
gauge boson, it will no longer annihilate into gauge boson pairs. Annihilation into matter 
fields from s-channel exchange of heavy gauge bosons is suppressed by the larger gauge boson 
mass (and depends on the number of channels into which the gauge boson can decay). This 
can allow for thermal dark matter that is somewhat lighter than in the doublet case, which 
can provide a better IDM fit to the direct detection data. 

4.2.2 A Light C/(l)^ Model 

IDM scattering at DAMA can also be mediated by a new light gauge boson that couples 
directly to the DM, but only very weakly to the visible SM sector. This can arise from 
an exotic U{l)x gauge boson with a mass on the order of a few GeV that has a kinetic 
mixing with electromagnetism on the order of e ~ 10~^ — 10~^. We describe a minimal 
super symmetric realization of such a scenario in this section. However, in this scenario and 
others like it, new mass scales must be put in by hand. Thus, we find this scenario less 
compelling from the point of view of naturalness (not technical naturalness) but permissible 
so we include it as a logical possibility. 

Supersymmetry affords a natural setting for light exotic gauge bosons in sectors that are 
somewhat shielded from the source of supersymmetry breaking. The smaller supersymmetry 
breaking soft terms in the hidden sector can then induce symmetry breaking in that sector 
at a scale that is parametrically smaller than the electrowcak scale. However, the minimal 
model we describe below also contains supersymmetric mass scales whose origin requires 
further explanation. 

Our IDM model contains two pairs of SM-singlet chiral superficlds, a and along with 
H and H'^, such that each pair has vector-like charges under a new U{l)x gauge symmetry. 
The model also contains a pure singlet, S, uncharged under both the SM and U{l)x. Both a 
and a'^ as well as S are assumed to be odd under an exact unbroken Z2 symmetry. We take 
the superpotential to be 

W D n' HH^ + Maaa^+^MsS^ + XiSa^'H + X2SaH^, (4.43) 

where we assume ~ Mg ~ TeV and \j! ~ GeV. We assume further that the fields in 
this sector are shielded from supersymmetry breaking relative to the MSSM sector. This 
can hold if supersymmetry breaking is mediated to the visible sector by gauge mediation 
through messengers charged only under the MSSM gauge group, for example. 

In the context of gauge mediation, kinetic mixing with hypcrcharge then induces effective 
charges for the gauge messengers under the U{X)x symmetry that are suppressed by the 
mixing parameter. The resulting soft scalar masses in the t/(l)a; sector are thus on the order 
of ~ {Qx/ g'Y i^'ec-i where e ~ lO"^ — 10~^ is the kinetic mixing parameter and m\c is 
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the soft scalar mass of the right-handed selectron [981 EHl EHl HOOj - Upon running down to 
lower energies, the scalar soft masses for H and can be induced to run negative by way 
of large Yukawa couplings, generating VEVs for these fields on the order of a GeV [59l llOOj . 
The H and scalars can also be destabilized at the origin by the contribution to the U{l)x 
D-term potential from the MSSM Higgs fields, which obtain small U{l)x charges from gauge 
kinetic mixing. The hidden sector VEVs generated in this way will be on the order of ^/tv. 
If the scale of gauge mediation is relatively high, the U{l)x sector can also receive additional 
small soft breaking contributions from gravity mediation [T5l llOlj . 

Integrating out S generates the effective superpotential 

W^e/z Df^'HH'^ + M^aa'^-^ {a^Hf - ^ {aH^f - ^{a^H){aH% {AAA) 



From this we obtain the fermion mass splitting 

v2 , x9 lTTr\2 



(4.45) 



where in the second equality above we have assumed Ai = A2 = A and {H) = {H^). With a 
small amount of suppression from the couplings, this is of the right size for IDM. From the 
F terms of H and we also get scalar splittings 

Vf D (4-46) 

-^Ma [H'^f « - ^ Ma {Hf a'^ a* + h.c. 

The first two terms here are subleading, while the second two generate scalar mass splittings 
of the right size with a small amount of additional suppression from the couplings Ai,2- Note 
that either the fermion or the scalar can be the dark matter state, whichever is lighter, 
depending on the soft terms. 

While the superpotential of Eq. (14.431) is technically natural, this model does not give 
an explanation for why Ma and Mg are so much larger than /x'. A value of /i' on the order 
of a GeV can perhaps arise naturally from an NMSSM-like extension of the Higgs sector in 
this model, as in Refs. [59 | llOOj . It could also arise from a Giudice-Masiero [I02j coupling to 
supergravity in high-scale gauge mediation |101j . The larger masses Ma and Ms could also 
potentially be generated by an NMSSM-like extension coupling to a*^^^ and S |100j . although 
this possibility would be more complicated. 

The model described above is similar to the one presented in Refs. p3 | llOlj . There, the 
larger mass is also put in by hand, the smaller mass is related to the breakdown of a non- 
Abelian gauge symmetry near a GeV, and the inelastic mass splitting arises radiatively from 
gauge boson loops. However, because this mechanism for inelastic splitting works only with 
a non-Abelian gauge group, the construction also requires a higher-dimensional operator in 
order to generate kinetic mixing with hypercharge, Tr((9*^"')" 14^°^,) B^'^/M'^, where O^"^^"" is 
a chiral adjoint operator under the exotic gauge group of dimension n. It requires additional 
model ingredients to explain the origin of the scale M which cannot be too large in order to 
obtain a large-enough mixing angle. 



36 



5 Conclusions 



In this paper we have considered the possibihty that DAMA is a true discovery of dark 
matter, and investigated the properties a theory of dark matter needs to have in order to 
account for the data. We have shown that inelastic dark matter is consistent with the findings 
of DAMA as well as other direct detection experiments and seems to provide a better account 
for the sum of this data than other proposed explanations such as light elastic dark matter 
or dark matter scattering off detector electrons. Extending the study of Ref. [Ilj, we find 



that heavier inelastic dark matter can give a reasonable fit to the data, particularly for lower 
values of the galactic DM escape velocity. 

An intriguing additional observation is that if the inelastic dark matter candidate is an 
electroweak doublet, it can simultaneously have the correct thermal relic abundance and a 
nucleon scattering cross section mediated by the SM in the range consistent with DAMA 
and other experiments. For a scalar doublet, this occurs when its mass is close to 525 GeV, 
while for a fermion doublet the mass should be about 1080 GeV. This makes electroweak- 
doublet inelastic dark matter candidates particularly attractive. 

The DAMA signal can also be explained by the inelastic scattering of DM off nuclei 
mediated by a new massive U{l)x gauge boson. This exotic gauge boson can either be heavy 
and couple directly to the SM, or very light and hidden. In the heavy case, the new gauge 
boson must be somewhat leptophobic and have large couplings to quarks and the DM. Light 
gauge bosons can work if they have a GeV-range mass and couple to the SM through a 
small kinetic mixing with the photon. In contrast to the case of electroweak doublet DM, 
however, the correct nucleon scattering cross section and thermal relic density do not arise 
automatically, and must be arranged by hand. 

Given that inelastic dark matter gives a compelling explanation for the DAMA result, 
it is of interest to understand what kind of models of IDM might work. The properties we 
need are clear. We need a particle of mass ~ 100 — 1000 GeV whose real and imaginary (or 
Weyl) components are split by about ~ 100 keV. Such a small splitting violates a U{1)dm 
global symmetry that is preserved by a Dirac or complex scalar mass term. This makes such 
candidates technically natural. 

We then address the question of overall naturalness. That is, why should the mass 
splitting be six orders of magnitude smaller than the overall mass scale? We found several 
candidate models. The first two work in the context of warped extra dimensions, where 
U{1)dm symmetry breaking can be sequestered. The symmetry breaking resides either in 
a large fermion Majorana mass or a small scalar holomorphic mass for a bulk singlet, that 
then mixes with an SU (2) l doublet dark matter candidate. The downside of these models 
is that we do not yet know if a warped extra dimension exists, or if the low UV scale needed 
in the scalar model is present. 

The remaining models work in the context of supersymmetry. In the most compelling 
supersymmetric model, the smallness of the splitting arises through the mixing between 
an SU{2)l doublet fermion dark matter candidate and a singlet which directly couples to 
U{1)dm breaking. The smallness of the mixing is attributable in part to large tan/3. The 
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downside is that we do not know if weak-scale supersymmetry is present, or why there 
should be a minor conspiracy of small numbers to give the necessary suppression. Other 
super symmetric models can work, but they generally require more assumptions or more 
complications to sufficiently isolate U{1)dm symmetry breaking. 

We do not yet know if the DAMA signal will prove to be new physics. If it is, we 
conclude that inelastic dark matter particles are excellent candidates. Reasonable (but not 
completely obvious) assumptions then lead to acceptable models of inelastic dark matter. 
We are very fortunate in that upcoming results from XENON, CRESST, and other direct 
detection experiments should be able to help determine if such models are likely correct. 
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A Appendix: Decoupling of KK Modes in the Warped 
Fermion Model 



In this appendix we show that the effect of the bulk singlet KK modes can be neglected 
when calculating the DM mass splitting in the warped fermion model of section 14.1.11 Note 
that this is despite the fact that the zero mode and KK modes have a large mixing due to 
the Majorana mass terms {i.e., Aqq ~ Amn)- The key point is that the KK modes all have 
vector-like U{1)dm preserving masses, and are very inefficient at communicating U{1)dm 
breaking. 

One can see this explicitly in the case of the first KK mode by including the couplings to 
the dark matter doublet as part of the mass matrix, and solving for the eigenvalues in the 
limit of a small Higgs VEV. For simplicity we can also neglect the Dirac mass of the doublet, 
which plays no role in the communication of U{1)dm breaking. The mass matrix is then 



1 



Qlc 



( ^00 

^01 





Am 
An 
mi 





m\ 





Civ 


/ 



+ h.c. 



(A.l) 



and the eigenvalues are determined from the roots of the characteristic polynomial Det(M 
XI) 

= X^-(Aoo + Au)\'-{Cy + Cy + mj)\^ 

+ (AooK + + AiiCy - 2AoiCoCiv'') A + mlcy, 



(A.2) 
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where we have used the fact that Aqo^h — ^oi — 0- hmit that f ^ goes to 0, there 

is a zero eigenvalue corresponding to the doublet. Thus we expect that the eigenvalue is 
proportional to v"^. Plugging in A = xv"^ and dropping terms of 0{v'^), we obtain 

^ xAoomjv^ + m\Clv^ + 0{v^) (A.3) 

so a; ~ and we see that at leading order the Majorana mass of the doublet is A ~ . 
This only depends on Aqo, and is independent of An and ^oi as we assumed in the estimates 
of section I4.1.1I 

The leading corrections to this formula are proportional to ~ and are suppressed 

2 

compared to the contribution from the zero mode by a factor ~ We have numerically 
checked that the effects of including more KK modes are also similarly suppressed, and that 
the sequence rapidly converges after the first few modes. We also note that this behavior is 
consistent with the results of Ref. I83l . 



B Appendix: Decoupling of Higher KK Modes in the 
Warped Scalar Model 

In this appendix we will show that the effective coupling of singlet KK modes to the DM 
doublet C^^^ becomes suppressed as ~ l/n for large n in the warped scalar model of 
section I4.1.2[ 

Substituting the form of the KK mode wavefunctions Eq. fl4.2ip into Eq. fl4.35p . we see 
that computing C^^^ requires performing the integral 



A 



A 



Close to y = ttR, the integrand is completely dominated by Jq. For large argument (large 
n), this Bessel function can be approximated as 



Changing variables to 2; = ^e^^ then gives 

rz+ 3 / 71 1 \ 

Cn" - kA J dz z-^ COS imnz - -{a + -) \ , (B.3) 
where z± = leH'^R^^)_ The cosine rapidly oscillates over the brane thickness at large n, 

3 

while the z~2 piece is relatively stable. Approximating the stable piece by its central value 



and performing the integral, one obtains 



■'n 



X^/ke 



-2iTkR 



TT 1 \ / TT 1, 

sin I m„2;i [a H — ) — sm m„2;_ [a + 



2' 2') \ 2' 2- 
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(B.4) 



and then one can place an approximate upper bound the magnitude of C^-^-^ , 

IQ^^I < . • (B.5) 

Thus, we see that the magnitude of C^^ falls roughly as ~ 1/n as we set out to show. 
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